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Abstract 

(N 

| Unintegrated gluon densities in nuclei, dipole-nucleus cross-sections and quark densi- 

ties are numerically investigated in the high-colour limit, with the scattering on a heavy 
nucleus exactly described by the sum of fan diagrams of BFKL pomerons. The initial 
condition for the evolution in rapidity is quickly forgotten, and the gluon density presents 
a " supersaturation" pattern, as previous studies indicated. Both dipole-nucleus cross- 
sections and quark densities present the expected saturation features. Identifying the 
position in transverse momentum I of the maximum of the gluon distribution with the 
saturation momentum Q s (y,b), at large rapidities all distributions depend only on one 
Ph| variable [l/Q s (y,b)] or [rQ s (y,b)]. 

Oh! 

J|h ; 1 Introduction 

In view of the current and forthcoming experimental investigations of the strong interaction 
with heavy nuclei at high energies, much attention has lately been devoted to the theoretical 
study of parton distributions inside a heavy nucleus at small values of the scaling variable x. 
It turned out that a particularly transparent approach follows from the colour dipole picture 
[1,2], in which the interaction of a probe with a target is presented via the interaction of the 
latter with a colour dipole, convoluted with the distribution of colour dipoles in the probe. 
In such a picture the fundamental quantity is the cross-section cr(Y, r) for the interaction 
of a colour dipole of a given transverse radius r with a target at a given rapidity Y. Much 
popularity has obtained the idea of " saturation" , which in terms of the cross-section a implies 
that at high energies the cross-section o~(Y, r) tends to a constant independent of r [3]. A 
particular ansatz chosen in [3] for the scattering of a dipole on the proton, 



with do — 23 mb and Rq diminishing with increasing rapidity Y, leads to a good description 
of the DIS data for the proton below x = 0.01. 

A less phenomenological treatment can be applied for the heavy nucleus target. In the 
framework of the colour dipole model it follows that, in the high-colour limit N c — > oo, the 
scattering on a heavy nucleus is exactly described by the sum of fan diagrams constructed of 
BFKL pomerons, each of them splitting into two. The resulting equation for the colour dipole 
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cross-section on the nucleus [4-7] was numerically solved in [7]. The gluon density introduced 
in [7] revealed a " supersaturation" behaviour, tending to zero at any fixed momentum k as 
Y — > oo. As a function of In k it proved to have a form of a soliton wave moving to the right 
with a constant velocity as Y increases. A more ambitious project is currently developed by 
the McLerran group, which can, in principle, lead to a description which does not employ the 
large number of colors limit [8-10]. Admitting that such an improvement is highly desirable, 
we think that the N c — > oo approach is much more feasible and can give a clear hint on the 
qualitative behaviour of parton densities and cross-sections for the heavy nucleus target. 

In this paper we continue studying the numerical solution of the BFKL fan diagram 
equation started in [7], with more precision and more attention to rapidities available at 
present or in the near future. We compare solutions obtained from a purely theoretical initial 
function (as in [7]) and from a phenomenologically supported one (like Eq. (1)). We find that 
the initial form is very quickly forgotten by the equation, so that at rapidities of the order 10 
the solution becomes independent of the chosen initial form. The behaviour of the solution 
at large energies is completely determined by the scale Q s (Y,b), depending on the energy 
and impact parameter b, at which the gluon density reaches its maximum value, so that both 
the parton densities and the dipole cross-section become universal functions of momentum 
or coordinate scaled with Q s . One may consider Q s as a "saturation momentum" introduced 
in [4,8,11]. The Y and b dependence of Q which follow from our numerical studies can be 
fitted by a simple formula: 

In Q s = a + cy + d \n[AT(b)} , y = ^±Y, (2) 

7T 

with 

c = 2.06, d = 0.62 ±0.07, 

Q s in GeV/c, and T(b) the nuclear profile function in (GeV/c) -2 defined so that / d 2 bT(b) = 
1; thus, AT(b) oc A 1 / 3 . Note that the value of c results lower than in our previous run devoted 
mostly to asymptotic energies and agrees with the predictions in [12] based on asymptotic 
estimates. The somewhat unexpected value of d seems to slightly depend on the choice of 
the initial function (hence the ±0.07). It implies that Q s oc A a , a ~ 2/9. The dipole cross- 
section on the nucleus at fixed impact parameter exhibits the expected behaviour, tending 
to unity at each r as Y — > oo. The quark density follows the pattern of saturation, tending 
to a constant value at small momenta, in full agreement with the predictions of [13]. 



2 The evolution equation and initial conditions 

As mentioned, in the colour dipole approach the cross-section of a probe (P) on the nucleus 
(A) is presented via the dipole cross-section: 

a PA (Y) = J d 2 rp P (r)a A (Y,r). (3) 

In its turn, the dipole cross-section on the nucleus is an integral over the impact parameter: 

a A (Y,r) = 2j d 2 b<Z>(Y,r,b), (4) 

where evidently 2$ has a meaning of a cross-section at fixed impact parameter. The evolution 
equation in Y can be most conveniently written for the function 

4>{Y,r,b) = -L$(y,r,6) (5) 
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in momentum space, where it reads [7] 

^- + Hbfkl^ <f>(y,q,b) = -0 2 (y,q,b), (6) 

with a rescaled rapidity y given in Eq. (2) and Hbfkl the BFKL Hamiltonian. 

Special attention has to be devoted to the initial function (f>(yo, q, b) = 4>o at y = yo, from 
which value one is starting the evolution. In the framework of the pure BFKL approach, 
with N c — > oo and A fixed, one should choose yo = (or any finite y) and take for (fio the 
contribution of the pure two gluon exchange with a single nucleon inside the nucleus. If A is 
large (A 1 / 3 of the same order as iV c ) then one should add all multiple interactions inside the 
nucleus, which sum into a Glauber cross-section: 

$ (r,6) = i_ e -^(6Vp(o,r) i (7) 

Here a p (0,r) is the dipole-proton cross-section generated by the two-gluon exchange: 

<r P (0,r) = \g* j d 2 r'G(0,r,r')p p (r'), (8) 

where p p (r) is the colour dipole density in the proton and G(0, r, r') is the BFKL Green 
function at Y = 0: 

rr' r < f r > 



G(0,r,r') = 1 + ln— , r >(<) = max(min){r, /}. (9) 

8ir r> V r<J 

Obviously the density p p is non-perturbative and not known. As in [7], to simplify the 
calculations we choose p p to be a normalized Yukawa distribution 

P P (r) = £ , (10) 

2ir r 

with p = 0.3 GeV adjusted to the nucleon radius value. With Eq. (10) we find 

AT(b)a p (0,r) = B [2C - 1 + 2 lnf - Ei(-f) (2 + f 2 ) +e~ f (l-f)] , r = pr, (11) 

where C is the Euler constant and dimensionless B carries all the information about the 
nucleus. For Pb at the center (b = 0) it reaches the value 0.12. Of course, this choice 
of p p may look rather arbitrary. We noted in [7] that calculations demonstrated a certain 
indifference of the evolution equation to the choice of the initial function at high enough y. To 
see this more clearly, in this run we also used an alternative initial function, more adjusted to 
the existing experimental data at comparatively small rapidities. A natural choice would be to 
take directly (1). However its analytic form makes it rather difficult to pass to the momentum 
space employed in our method of evolution. Therefore we choose a slightly different form for 
the phenomenologically motivated initial dipole cross-section on the nucleon: 

AT{b)af = B (l - e- f f , (12) 

with the same parameters as in [3] and Yq corresponding to x = 0.01. The cross-section in Eq. 
(12) has the same asymptotic behaviour as that in Eq. (1), both at r — > and r — > 00. For 
finite r it is slightly different, but this difference is not too large and, after all, the choice (1) is 
also quite ad hoc. Note that the asymptotic behaviour of the phenomenological cross-section 
is different from that of Eq. (11), which both at r — > and r — > 00 contains an extra | lnr| 
factor. The dimensionless B in Eq. (12) results substantially larger than for Eq. (11) and is 
of the order of 5 for Pb at the center. 

However, as we shall presently see, in spite of the difference both in the asymptotic 
behaviour and in the overall normalization, starting from y ~ 2 the results of the evolution 
begin to look practically identical for both choices of the initial function, except for the overall 
scale Q s mentioned in the Introduction. 
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3 The gluon density and dipole cross-sections 



We define the gluon density as in [7] : 

d[xG(x,k 2 ,b)} 27V 2N 

fflwp = ^2 k V kHv, k, b) = ^Hy, k, b), (13) 

with y = (N c a s /Tr) ln(l/x). This definition follows the logic of the BFKL approach in which 
it corresponds to the average of two gluon fields in the nucleus, in the axial gauge adopted 
in this approach. It also naturally appears in the expression for the structure function of the 
nucleus (see [7]). There exist different definitions of the gluon density in the literature (see 
[13]). We are not going to discuss the problem of a "correct" gluon density here, since in any 
case it is not a directly measurable quantity, but rather serves to calculate the latter (see 
comments in [14]). From a pragmatic point of view our definition is well supported, since all 
observables can be directly related to (13). 

The results of our evolution for the gluon density are shown in Figures 1 and 2. Not to 
bind ourselves to a particular value of the coupling constant a s , we rather present h(y, k, b) 
as a function of the rescaled rapidity y. The actual gluon density at physical rapidity is 
obtained after rescaling both h and y according to Eqs. (2) and (13). 

In Fig. 1 we show the gluon densities at the first stage of the evolution up to y = 
1, starting from the initial functions corresponding to the Glauberized two-gluon exchange 
contribution (11) (theoretical initial function, TIF) and to the phenomenological cross-section 
(12) (phenomenological initial function, PIF), respectively. One observes a large difference 
between the two at the beginning of the evolution, which, however, is gradually disappearing. 
Starting from y ~ 2 the form of the gluon distribution becomes practically identical for both 
initial functions. We introduce Q s (y,b) as the the momentum at which the density reaches 
its maximum. For the two initial functions a numerical fit gives Eq. (2) with a = —2.21 
(-0.88) , c = 2.09 (2.04) and d = 0.69 (0.55) for TIF (PIF). Plotted against k/Q s (y,b) the 
gluon distributions for both initial functions and for all values of b and y > 2 practically fall 
onto a universal curve. The degree of universality is illustrated in Fig. 2 where we present 
the gluon densities for both TIF and PIF at y = 2.4 and 4, and at b = and b = 0.98Ra, for 
the Pb target. Some differences in the curves are certainly visible, which however diminish 
as y grows. The form of the universal curve can be fitted by an expression 

h(0 = 0.295 exp "ggg ), £ = lnfc, £o = mQ s (y,6) (14) 

(all momenta in GeV/c), from which one observes that the density falls as a function of 
momentum both at k — > and k — ► oo. 

The gluon density (13) is trivially related to the dipole cross-section on the nucleus at 
fixed b: 

*{y,r,b) = l- I ^h(y,k,r)^. (15) 



2irk 2 

So, to determine <&, all one has to do is to perform a Bessel transform of h/k. Fig. 3 shows 
the dipole cross-sections at low rapidities. At larger rapidities y > 2 the scaling properties 
of h imply that <I> becomes a universal function of rQ s (y, b). Its form is illustrated in Fig. 4, 
where we show <3>(r) for Pb at b = and y = 3.4 (for PIF). It can be fitted with a formula 



<j>( r ) 



1 — exp ^— uj 2 r 2 



i/S 



uj = 198.5 GeV/c, 5 = 5.48. (16) 
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4 The quark density 



The definition of the quark density can be taken from [13]: 

d[xq(x,l,b)\ _ a s Q 2 



d 2 ld 2 b 
,bib 2 



(2tt): 



■ f 1 dad^cfhe-^-^ 
Jo 



(a 2 + (1 - a) 2 ) e 2 -L-lK 1 (e6 1 )K 1 (e6 2 ) + 4QV(1 - a) 2 K (e&i)K (e& 2 ) 
d 2 k 1 d[xG{x,k,b)} 



I 



(2vr) 2 k 2 d 2 bdk 2 



1 + e 



-ik(bi-b 2 ) 



-ikbi 



_ikb2 



(17) 



with e 2 = Q 2 a(l — a). This definition is based on the form of the interaction with the target 
of a virtual current which splits into a qq pair. For the small x region one may also raise 
objections as to its physical meaning, since then the interference diagram gives a non-zero 
contribution [14]. However, for lack of a better definition, we shall use Eq. (17). 

Performing part of the integrations and using Eq. (13), we express the quark density as 



d[xq(x, I, b)] 
d 2 W 2 b 



= 7r 4 f(y,l,b), 



(18) 



where y = (N c a s /ir)\n(l/x) and function /, independent of the coupling constant, is given 
by an integral over momenta: 



f 

Jo 



dk 

T 



h(y,k,b)I(l,k). 



Here I is a sum of transversal (T) and longitudinal (L) parts, with 



I T 



Q 2 C 
Jo 



da 



a 2 + (l 



a) 



2 + k 2 



( e 2 + / 2 + £; 2 ) X 3 -e 2 A 2 



and 



where 



I L = 4Q 4 [ da 2 a{l - a) 2 [(e 2 + I 2 + A: 2 ) X 3 + A 2 - 2 X X 



e 2 + I 2 ' X ^(e' 2 + (l + ky 2 )(e 2 + (l-k)' 2 ) 



(19) 



(20) 



(21) 



(22) 



Note that the quark density defined by Eqs. (17)-(22) depends (weakly) on the virtuality 
of the probe Q 2 . Physically meaningful results correspond to the limiting case Q 2 —> oo. 
In this limit the longitudinal part of / gives no contribution and in the transversal part the 
integration over a is reduced to that over e 2 : 



Q 2 C 
Jo 



da 



a 2 + (l-af 



Jo 



de 2 , 



(23) 



so that / and / become independent of Q 2 . 

These formulas give the quark density in terms of the gluon density, that is, in terms of 
the function h, which we have found numerically. Doing the integration over k in Eq. (19) 
also numerically, we obtain function / and therefore the quark density related to it by a factor 
(see Eq. (18)). Our results for the quark density (in fact for /) are presented in Figs. 5 and 
6. As before, in Fig. 5 we show the quark densities at the beginning of the evolution for the 
initial functions TIF and PIF respectively. At higher y > 2 and due to scaling properties of 
the gluon density, the quark density also becomes a universal function of l/Q s (y, b). Its form 
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is illustrated in Fig. 6, where we show /(/) for Pb at b = and y = 3.6 (for PIF). It can be 
fitted by 



1 - exp 



I - 2 



f(t) = —Z-± ^ J -L- t = 2[ , 3 = 2.44. (24) 

J W 2vr 8 t ' V 511.0 GeV/c/ ' H K ' 

Note that at I — > the quark densities tend to a constant value in agreement with the 
prediction of A. H. Mueller [13]: 



d[xq(x, I, ft)] N c 

d 2 id 2 b ^° 2^4 ' 

which implies 

f(l) — 1.58 • 10~ 4 ; 
at I — ► oo they acquire a perturbative form (oc l/l 2 ). 



(25) 



5 Discussion 



The new run of numerical investigations of the BFKL fan diagram equation for the gluon 
density in heavy nuclei fully confirms our previous conclusions made in [7]. The key one is 
that at rapidities of the order 10 the density forgets its initial form and becomes a universal 
function of momentum scaled by the saturation momentum Q s (y,b). The latter grows as a 
power of energy and as A a with a of the order of 2/9. Both the dipole-nucleus cross-section 
and the quark density behave in agreement with a general idea of saturation, tending to a 
known constant (unity for the dipole cross-section at fixed impact parameter) as rapidity 
grows. 

Numerical studies, as experimental ones, need confirmation from independent groups. 
To our knowledge, up to now there appeared only two papers devoted to the solution of 
the non-linear BFKL equation, both of them for the nucleon target. In [15] a simple Pade 
technique is used to solve the non- linear equation obtained in [6]. Employing for the linear 
case a solution which embodies both BFKL and DGLAP behaviour and a running coupling 
constant, the authors present results for the integrated gluon distribution at Large Hadron 
Collider energies, finding a suppression effect, due to the non-linearity, of a factor ~ 2. More 
appropriate for comparison with our results are those of [16] , in which an iteration technique 
in coordinate space was used and results for the dipole cross-sections were shown. The chosen 
initial function is different from both our choices. It attempts to include the DGLAP evolution 
at the initial stage. We are not going to discuss here the viability of such an attempt for 
the nucleus target. Still our results show that at rapidities studied in [16], of the order of 
12-r- 15, which with a s = 0.25 mean y ~ 3.3 4- 3.8, the dipole cross-sections <3?(r) (N(r) in the 
notation of [16]) should become a universal function of f = rQ s . The comparison indicates 
that although the general behaviour of $ found in [16] well agrees with ours, the universality 
in the f-behaviour is not observed there. The slope of their curves in In f is steeper than for 
our universal curve and this difference grows with rapidity^]. 

1 After the first version of our work, an extension of the one in [16] to the nuclear case has appeared, [17]. 
As in [16], the initial conditions are different from ours, but the authors find a dependence Q s oc A a with 
a — 1/6 at very small x, quite similar to the one we obtain. Besides, they claim that the universality of gluon 
distributions and dipole cross-sections is also observed in their numerical solutions. 
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Figure captions 

Fig. 1: The gluon densities at the first stage of the evolution y < 1, for Pb target at b = 0. 
Solid and dashed curves show the densities evolved from TIF and PIF at y = respectively. 
Curves from left to right correspond to y = 0.0, 0.4 and 1.0. 

Fig. 2: Gluon densities at y = 2.2 and 4.0, for Pb target at b = and b = 0.98Ra, both for 

TIF (solid curves) and PIF (dashed curves) as initial functions, plotted against k/Q s (y,b). 

Fig. 3: The dipole cross-sections <J>(r) at the first stage of the evolution y < 1, for Pb target 

at b = 0. Solid and dashed curves show the cross-sections for TIF and PIF as initial functions 

respectively. Curves from right to left correspond to y = 0.0, 0.4 and 1.0. 

Fig. 4: The dipole cross-sections <3?(r) at y = 3.4 for Pb at b = 0, with PIF as initial 

function. 

Fig. 5: The quark densities f(l) at the first stage of the evolution y < 1, for Pb target at 
6 = 0. Solid and dashed curves show the quark densities for TIF and PIF as initial functions 
respectively. Curves from left to right correspond to y = 0.0, 0.4 and 1.0. 
Fig. 6: The quark densities f(l) at y = 3.6 for Pb at b = 0, with PIF as initial function. 
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Figure 3: 
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